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Bacterial attachment to a fluid interface can lead to the formation of a film with physicochemical 
properties that evolve with time. We study the time evolution of interface (micro)mechanics for 
interfaces between oil and bacterial suspensions by following the motion of colloidal probes trapped 
by capillarity to determine the interface microrheology. Initially, active bacteria at and near the 
interface drive superdiffusive motion of the colloidal probes. Over timescales of minutes, the bacteria 
form a viscoelastic film which we discuss as a quasi-two-dimensional, active, glassy system. To study 
late stage mechanics of the film, we use pendant drop elastometry. The films, grown over tens of 
hours on oil drops, are expanded and compressed by changing the drop volume. For small strains, by 
modeling the films as 2D Hookean solids, we estimate the film elastic moduli, finding values similar 
to those reported in the literature for the bacteria themselves. For large strains, the films are highly 
hysteretic. Finally, from wrinkles formed on highly compressed drops, we estimate film bending 
energies. The dramatic restructuring of the interface by such robust Aims has broad implications, 
e.g. in the study of active colloids, in understanding the community dynamics of bacteria, and in 
applied settings including bioremediation. 


I. INTRODUCTION 

There is intense interest in the dynamics of microbes in 
suspensions, as they can be considered as active colloids. 
Bacteria in planktonic states convert chemical energy to 
propel themselves, and in dense suspensions, display rich 
collective behavior [1-3]. In this context, swimming mi¬ 
crobes have inspired studies of reductionist systems in 
which particles convert chemical fuel into motion [4, 5], 
which are analyzed in terms of non-equilibrium statisti¬ 
cal mechanics [6]. However, microbes are not merely col¬ 
loids that move; they are living entities with complex sur¬ 
face chemistries and structures, which allow them to form 
complicated networks and to respond to external stimuli. 
Microbes also produce and secrete a variety of materi¬ 
als ranging from small molecule surfactants, polysaccha¬ 
rides, and growth factors that allow them to signal each 
other and to exist in communities [7]. Via these secre¬ 
tions and associations, bacteria restructure their physic¬ 
ochemical environment. 

The most common form of communal behavior is biofilm 
formation; an estimated 99% of the world’s bacteria 
are sequestered in biofilms in various stages of growth 
[8, 9]. Biofilms comprise adherent cells in a com¬ 
plex three-dimensional matrix of extracellular polymeric 
substances[10] which resemble a disordered colloidal sus¬ 
pension in a polymer solution or gel [11]. These complex 
multicellular 3D macrostructures on solid surfaces grow 
in stages; microbes adhere to a surface, and grow into 
microcolonies that eventually form thick Aims [12]. To 
form biofilms, microbes exploit synergistic interactions 
between cells coordinated by signaling via secreted en¬ 
zymes and proteins [13]. In nature, the myriad roles 
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played by biofilms continue to be identified. For exam¬ 
ple, biofilm formation can be beneficial for community 
survival, as biofilms protect the growing population of 
microbes against potential physicochemical and biologi¬ 
cal attacks [14, 15]. Biofilms can protect bacteria from 
antimicrobial agents, with implications in antibiotic resis¬ 
tance [9, 16], and can confer evolutionary advantages via 
cooperation of co-colonized bacteria strains [17]. These 
protective mechanisms can be compromised, however, as 
motile bacteria of competing species can exploit the per¬ 
meable, open structure of biohlms to invade and displace 
an established biofilm colony [18]. 

In more applied settings, biohlms play important roles. 
They can be remarkably strong elastic structures, and 
grow readily on bounding surfaces of pipes, reactors, 
ships, etc. As such, they decrease the efflux from fac¬ 
tories, foul Alters and bioreactors [19], and are present 
in drag-increasing biocolonies on ships [20]. Biohlms can 
also be beneAcial; for example, in bioremediation, bacte¬ 
ria and their biohlms remove toxins from their surround¬ 
ings [21, 22]. 

The formation, structure, and mechanics of biohlms 
supported on solid surfaces are relatively well stud¬ 
ied [11, 23-25]. In contrast, Alms formed by microbes 
trapped at fluids interfaces have been less explored. We 
refer to such Aims involving bacterial cells as /ilms of 
bacteria at interfaces (FBI). The nature of FBI and 
their relationship to the better-studied solid supported 
biohlms remain to be elucidated in terms of FBI com¬ 
position, dynamical behavior and biological roles. There 
are a few notable recent studies in this held; for example, 
at air-aqueous interfaces, the structure and properties of 
interfacial biohlms or pellicles of B. subtilis [26, 27] have 
been studied to investigate the role of amyloid flbers in 
providing structural integrity to the Aims. The time evo¬ 
lution and mechanics of FBI formed by E. coli have been 
also studied, revealing rich interface rheology dependent 
on microbial secretions [28-32] with implications in infec- 
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tion [29, 32], FBI at oil-water interfaces have been stud¬ 
ied [30, 33-38], motivated by their relevance to petroleum 
technologies, including oil recovery from oil sands [33- 
38]. Such studies are also relevant to bioremediation, 
e.g. in the context of oil-consuming bacteria that appear 
near oil spill sites. FBI at oil-water interfaces form elas¬ 
tic skins, with film elasticities that depend on the surface 
structures of the microbes [33, 34, 39]. 

Here we undertake a comprehensive study of the dynamic 
and spatial dependence of FBI formation in terms of 
film (micro)mechanics. We study FBI formation at ini¬ 
tially bacteria-free hexadecane-aqueous interfaces for the 
model organism Pseudomonas sp. ATCC 27259, strain 
P62. Our experiments reveal three stages of behavior: 
active films, viscoelastic films and elastic films. Initially, 
bacteria near the interface can be highly motile, forming 
an active layer. Over time, a film of adherent bacteria 
trapped in the interface forms; in the absence of nutrient 
addition, these bacteria are typically non-motile. At this 
stage, the FBI is viscoelastic. Particle-tracking experi¬ 
ments, in which film properties are revealed by trajecto¬ 
ries of colloidal probes trapped at the interface, provide 
information about these stages. 

To our knowledge, this is the first study that captures 
all of these stages, including the observation of active 
bacteria trapped at fluid interfaces and their transition 
to a viscoelastic film. The behavior of confined suspen¬ 
sions of motile bacteria near solids has been intensely 
studied. Near solid boundaries, bacteria swim in recir¬ 
culating patterns influenced by their confinement, their 
shape, and the nature of their bounding surfaces[40, 41]. 
The emergence of these patterns is attributed to hydro- 
dynamic interactions with shear flows [42] and with image 
charges that enforce the no slip boundary condition[43- 
45]. We characterize the dynamics of our active layer by 
analyzing the motion of colloidal probes in the interface. 
Bacteria trapped at the interface by capillarity influence 
the motion of these probes. A population of non-adherent 
bacteria also influence the interface dynamics, leading to 
collective behavior which we compare to prior work in 
quasi-2D systems. 

We infer that the viscoelastic transition is related to 
the secretion of polysaccharides and surfactants in the 
interface. This work complements existing studies in 
which secreted surfactants have been reported to play 
conflicting roles, as secreted surfactants are implicated 
in immobilization[46], or in inducing superdiffusive mo¬ 
tion at the air-water interface [2]. 

We address the FBI micromechanical properties from a 
fundamental perspective in the context of active, soft- 
glassy systems. FBI eventually form elastic, solid films. 
We probe the dilatational and bending moduli of these 
films using imposed area perturbations via pendant drop 
elastometry [47, 48], and provide evidence of nonlinear 
and hysteretic behavior for large area strains. These re¬ 
markably robust films significantly modify the mechanics 
of the interfaces on which they form. This restructuring 
of the physicochemical environment has broad implica¬ 


tions in nature and in applied settings. 

II. EXPERIMENTAL METHODS 
A. Sample preparation 

Pseudomonas sp. ATCC 27259, strain P62, chosen 
as a model organism, are cultured in ATCC Medium: 
3 Nutrient Broth in 24 ppt Instant Ocean to simulate 
a middle-level marsh salinity. The cultures are grown 
on a table-top shaker at 150 rpm at room temperature 
for 24 h to mid-to-late exponential phase. For the FBI 
characterization experiments, a bacteria suspension free 
of surface active proteins from the nutrient broth is 
prepared with a typical washing protocol[34, 49] by cen¬ 
trifuging the culture for 5 minutes at 5000 x g, decanting 
the nutrient broth or supernatant and re-suspending the 
pellet in Instant Ocean three times. 


B. Particle tracking 

The particle-tracking measurements are performed in a 
1.0-cm I.D. cylindrical vessel with an inner surface whose 
bottom half is aluminium and top half is Teflon. When 
the cell is filled to the appropriate level, the aluminium- 
Teflon seam pins the oil-water interface, creating a flat 
interface with no meniscus. The base of the cylinder 
rests on an untreated glass coverslip, sealed with silica 
vacuum grease. To begin each experiment, the cylin¬ 
der is filled nearly to the seam with 0.5 ml of Instant 
Ocean. A spreading solution containing charge-stabilized 
polystyrene spheres (Invitrogen) with radius Rs = 0.5 
g,m in a mixture of equal parts by volume water and 
isopropanol is prepared in advance and sonicated to dis¬ 
perse any colloidal aggregates. To introduce the colloids 
to the interface, a 10 /xL droplet of the spreading solu¬ 
tion is gently placed in contact with the Instant Ocean 
surface. The solution wets the surface, and the colloids, 
which are slightly denser than water but hydrophobic, 
disperse across the interface. Droplets of hexadecane are 
then promptly placed on this system to form a film ap¬ 
proximately 2 mm thick which covers the colloids and 
the aqueous suspension. The age of the sample ta is 
measured from the instant of oil-aqueous interface for¬ 
mation. After each experiment, the cylindrical vessel is 
cleaned thoroughly by scrubbing and sonicating in Al- 
conox soap solution, acetone, and isopropanol, and then 
rinsed repeatedly in deionized water. Control microrheol¬ 
ogy experiments of dilute nutrient broth were performed 
absent bacteria to confirm the evolution of the interface 
was due to bacteria and film development. We observe 
the colloids at the interface using an upright bright-held 
microscope with a 50x objective. A camera (Zeiss Ax- 
iolmager Mlm) records a 200 gm by 150 /xm held of 
view at 60 frames per second, which set the shortest 
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lag time over which probe motion can be characterized. 
Probe trajectories are extracted from the video using a 
custom Python implementation [50] of the widely used 
Crocker-Grier multiple-particle-tracking algorithm [51]. 
Static and dynamic errors in the particle tracking are 
taken into account to avoid distortion in the measure¬ 
ments of the particle trajectories [52]. Trajectories are 
extracted from video segments short enough in duration 
(typically 1-2 minutes) that no apparent change in par¬ 
ticle mobility owing to the evolving interfacial properties 
can be discerned. Typically 30-200 probes are in view at 
a time, constituting up to 0.3% surface coverage. 


C. Pendant drop elastometry 


Pendant drop elastometry is used to study the evolu¬ 
tion of the tension of the oil-water interface as the FBI 
forms. Changes in tension are also recorded under con¬ 
trolled deformations to the film achieved by expanding 
or contracting the droplet. These stress-strain observa¬ 
tion are compared to appropriate models to extract ma¬ 
terial properties. A droplet of hexadecane is formed, in 
contact with the bacterial suspension, at the tip of an 
inverted needle by injection from a syringe. The drop is 
positioned in the path of a beam of light which projects 
the drop silhouette onto a CCD camera. The tension 
7 o is determined by comparing the edge of the droplet 
to a numerical solution of the Young-Laplace equation 
for an isotropic interface Ap = 7 o(ks + K 4 ,) where Ap is 
the pressure jump across drop interface, and Ks and 
are the meridional and parallel curvatures, respectively. 
This form for the Young-Laplace equation can be recast 
in dimensionless form in terms of a Bond number for the 
droplet Bo = , where Rq is the radius of curvature 

at the drop apex, Ap is the density difference between the 
drop and the external phase, and g is the gravitational 
acceleration constant. A comparison of numerical solu¬ 
tions to this equation and digitized experimental drop 
edges allows the tension 70 to be determined for a given 
drop image. 

This analysis is appropriate immediately after drop for¬ 
mation. In later stages, the film behaves as a thin elastic 
sheet at the interface. If this sheet is unstrained, the 
isotropic form of the Young-Laplace equation and the 
standard pendant drop tensiometry remain appropriate. 
If, however, this sheet is strained, the principal stretches 
As and are not isotropic, nor are the tensions, given 
by Ts in the meridional direction and in the parallel 
direction. In this case, the relevant form for the Young- 
Laplace equation is Ap = TsKs + Tc/,Kc/) with the associ¬ 
ated tangential stress balance — -I- = 0, 

where s is the arclength measured from the apex of the 
droplet, 4 ) is the turning angle and r is the radial lo¬ 
cation of the drop edge in cylindrical polar coordinates. 
The drop shape is then determined by the relevant con¬ 
stitutive equation that relates local tensions (ts and r^) 
to local stretches (As and A^). For an axisymmetric 


droplet in the limit of small deformations, the film can 
be described by a 2D Hookean model, for which the elas¬ 
tic response is dependent on the two material parame¬ 
ters, the two-dimensional elastic modulus K 2 D and the 
two-dimensional Poisson ratio V 2 d- To determine these 
parameters, the drop is first imaged in its isotropic un¬ 
strained state to determine 79 . Thereafter, the drop is 
deformed from that state and its silhouette is captured 
using a CCD camera. The edge of the strained drop is 
compared to solutions of the anisotropic Young-Laplace 
equation (see, c.f. Carvajal et al. [48] and Knoche et al. 
[47] for details). To do so, we write an objective func¬ 
tion for the difference between the experimental and nu¬ 
merical profiles and minimize it in a first order Newton- 
Raphson scheme with respect to the location of the axis 
of symmetry of the drop, the Bond number in terms of 
the unstrained tension, and the material parameters K 2 D 
and 1^20 ■ From this analysis, we determine the elastic 
modulus K 2 d- We find, however, that this method is 
insensitive to the two-dimensional Poisson ratio i' 2 D- Fi¬ 
nally, following Knoche et al. [47], we estimate the bend¬ 
ing modulus for the film from the wavelength of wrin¬ 
kles that form under strong compression according to 

T 

Eb ~ ie^ 2£2 where Ac is the wavelength of the wrin¬ 
kles at the onset of wrinkling and is the length of the 
wrinkles which appear on the neck of the drop. 


III. RESULTS AND DISCUSSIONS 

The interface between hexadecane and the aqueous 
bacteria suspension evolves with age, changing measur¬ 
ably over a timescale of minutes. We discuss three qual¬ 
itatively distinct stages of the FBI development during 
which it can be characterized as (i) active, (ii) viscoelas¬ 
tic, and finally, (iii) as a solid elastic film. The first two 
stages of development are investigated primarily through 
the particle-tracking experiments, while the third stage 
of biofilm mechanics is characterized using the pendant 
drop experiments. 


A. Particle tracking 

Tracer Motion at an Active Interface 

Immediately following the formation of a fresh interface 
between oil and the bacteria suspension, bacteria were 
sometimes observed to associate with the interface 
and remain motile once attached. Bacteria attach to 
the interface in both end-on and side-on orientations. 
Often, bacteria at the interface were accompanied by 
a population of especially mobile bacteria just beneath 
the interface. In these cases, the colloidal motion at 
the interface was strongly affected by hydrodynamic 
interactions with the swimming bacteria at and near 
the interface and by direct collisions with those at the 
interface. Typically, the concentration of swimming 
bacteria embodied a 10 % area fraction of the interface - 
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large enough to lead to collective motion such as swirling, 
which influenced colloidal motion. A substantial litera¬ 
ture has addressed the mobility of colloidal probes in the 
presence of motile bacteria and other microbial swim¬ 
mers both in bulk (3D) and in quasi-two-dimensional 
contexts [3, 53-61]. One motivation for studying the 
colloidal dynamics in suspensions of swimming microbes 
is their utility as model systems for investigating the 
non-equilibrium statistical mechanics of active complex 
fluids. In addition, the colloid motion provides insight 
into biomixing, the enhanced transport of nutrients 
and other material by active suspensions. Since film 
formation relies on the transport of polysaccharides 
and other constituents to and within the interface, such 
biomixing could be an important feature of the film 
development. Hence, we have characterized the colloidal 
motion during this active phase in some detail. 

Fig. 1(a) shows the colloids’ ensemble- 
average mean-squared displacement (Ar^(<)) = 

+1) — where the average is over 

particles i and time t' during this active stage of FBI 
formation, and ri(t') is the time-dependent position 
of the ith particle. For reference, the mean-squared 
displacement of colloids at the oil interface of pure 
Instant Ocean containing no bacteria is also shown. In 
the absence of bacteria, (Ar^(t)) varies linearly with 
lag time t indicating simple diffusion, (Ar^(t)) = dUgt, 
with a diffusion coefficient, Dq = 0.122 consistent 

with the viscosities of water and hexadecane. The 
mean-squared displacement of the colloids at the active 
interface similarly varies linearly with t at large lag times, 
but with an enhanced effective diffusion coefficient. At 
smaller lag times, the mean-squared displacement grows 
more rapidly than linearly, indicating superdiffusive 
motion. Such superdiffusive motion is a common feature 
of colloidal motion in microbial suspensions and signals 
temporal correlations in the forcing of the colloids due 
to hydrodynamic interactions with the swimmers. A 
simple model for these correlations ascribes to them 
a single characteristic correlation time t, so that 
the particle velocities have an exponentially decaying 
memory, (vi(t') • Vi(F -|- t))i,t' ~ exp(—t/r) [53]. Such 
velocity correlations lead directly to a mean-squared 
displacement of the form [61], 

(Ar2(t)) =4D[t + T(e-‘/"-l)]. (1) 

In the limit of short lag times, t ^ t, this form predicts 
ballistic motion, (Ar^(t)) = at large lag times it 

reduces to diffusive motion, (Ar^(t)) = ADt, with effec¬ 
tive diffusion coefficient D. The solid line in Fig. 1(a) 
is the result of a fit using this form, which describes 
the data accurately and gives D = 1.8 ± O.l'^’^Ys 
r = 0.053 ± 0.007s. The value of D in relation to the 
diffusivity in the absence of bacteria, ^/dq ~ 15, indi¬ 
cates that biomixing strongly influences the interface dur¬ 
ing this early stage of Him formation. We note further 
that this factor likely underestimates the enhancement 


in tracer mobility due to the swimming bacteria since, as 
described below, even at these early interface ages the in¬ 
cipient film imparts an interfacial viscosity that reduces 
thermal diffusivity. 

The success of Eq. (1) in capturing the form of (Ar^(f)) 
is consistent with several earlier studies of colloidal mo¬ 
tion in active microbial suspensions and hence indicates 
that the velocities of tracers in such suspensions indeed 
have exponentially decaying memory. However, we can 
interrogate these correlations more directly by defining 
an instantaneous direction of motion for each particle. 


ni(T) 


d(t) - D(r - 5t) 
|d(t) - ri(r - (5 t)1 ’ 


( 2 ) 


where 5t = i/go s is the time between successive video 
frames and i labels the particles, and by examining its 
time-time autocorrelation function. 


4>„(t) = (n,(T-kt) • n,(T)),_^ (3) 

The significance of <I’„(t) is that it quantifies the per¬ 
sistence in the direction of the colloids’ trajectories. As 
shown in Fig. 1(b), 4>„(t) is effectively zero at t > 0.15 s, 
setting the typical time required for the colloids’ direction 
of motion to randomize completely. Notably, 4>„(t) does 
not decay monotonically but has a peak near t = 0.03 s. 
We attribute this peak to a tendency for colloids to fol¬ 
low U-shaped trajectories on short time scales due to 
hydrodynamic interactions with bacteria swimming past 
in close proximity [59], and hence to make negative con¬ 
tributions to $„(t). At larger times, <I’„(t) decays ex¬ 
ponentially, as shown by the line in Fig. 1(b), which is 
an exponential fit to the data. The correlation time ob¬ 
tained from the fit is 0.075±0.01s, which is in reasonable 
agreement with the characteristic time r for velocity cor¬ 
relations implied by fitting (Ar^(t)) with Eq. 1. Thus, 
through 4>„(f) we observe directly the nature of the cor¬ 
related tracer dynamics in an active suspension inferred 
by the analysis of (Ar^(t)). 

While the colloids’ Brownian dynamics at large lag times 
with diffusivity D suggest that the suspension of swim¬ 
ming bacteria acts like a thermal bath with large effec¬ 
tive temperature, several previous studies have empha¬ 
sized that the statistical properties of the colloidal dis¬ 
placements differ from those expected for a system in 
thermal equilibrium [56, 59, 61]. These differences are 
apparent in the probability distribution function (PDE) 
for displacements at fixed lag time Pt{Ax), where Ax is 
the displacement along one direction. Eig. 2(a) shows 
Pt{Ax) at t = 0.0167 s, t = 0.15 s and t = 1.3 s, lag 
times spanning the superdiffusive and diffusive behav¬ 
ior in (Ar^(t)). The PDEs of particles undergoing ther¬ 
mal diffusion in equilibrium would be Gaussian. In all 
cases, the PDEs in the active, incipient film show clear 
deviations from Gaussian forms, with tails at large [Ax] 
that signal enhanced probability of large displacements. 
Qualitatively similar non-Gaussian PDEs have been ob¬ 
served previously among tracers in microbial suspensions 
and have been associated with advection-enhanced large 
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displacements due to hydrodynamic encounters between 
the colloids and swimmers [59, 61]. To compare the 
form and magnitude of the non-Gaussian contributions 
to Pt{Ax) at different lag times more closely, we plot in 
Fig. 2(b) the normalized PDFs with displacement nor¬ 
malized by the root mean-squared displacement. Re¬ 
markably, the normalized PDFs collapse onto a single 
lineshape, indicating that the distribution function main¬ 
tains a self-similar form with increasing lag time. Such 
self-similarity is generally unexpected and implies par¬ 
ticular attributes about the colloidal dynamics, includ¬ 
ing (i) that the fraction of colloids in the non-Gaussian 
population remains constant over the range of lag times 
probed, and (ii) that the Gaussian and non-Gaussian dis¬ 
placements grow as the same function of lag time, first 
superdiffusively at short lag times and then diffusively 
at longer lag times. Self-similar distributions with non- 
Gaussian tails were also observed among tracer displace¬ 
ments within bulk (3D) suspensions of the eukaryotic 
microorganism Chlamydomonas[59]. However, in that 
case, (Ar^(t)) displayed diffusive behaviour over the en¬ 
tire range of lag times probed. The collapse in Fig. 2(b) is 
particularly notable because the lag times span both the 
superdiffusive and diffusive regimes. In both the previous 
case of tracers among swimming Chlamydomonas and in 
our case of tracers in an incipient FBI, the non-Gaussian 
contributions to the probability distribution function fol¬ 
low a Laplace distribution, so that the total PDF can be 
described as the sum of two parts, 

Pt{Ax) = -^exp[- ( —)2] + ^exp(-|Ax|/0, 

{2Tia^y' 2 cr ^ 

(4) 

where / is the fraction of colloids in the non-Gaussian 
distribution, ^ is the characteristic length of their dis¬ 
placements, and cr is the width characterizing the dis¬ 
placements of the Gaussian distribution. The solid line 
in Fig. 2(b) is the result of a fit to the data using this 
form. The strong similarity between self-similar PDFs 
of tracer displacements in the incipient FBI and those in 
bulk suspensions of Chlamydomonas is surprising since 
the tracer displacements in quasi-2D films of the Chlamy¬ 
domonas suspensions display qualitatively different non- 
Gaussian distributions, and the authors of those stud¬ 
ies attributed the difference to the differing fluid veloc¬ 
ity fields generated by the force dipoles of the Chlamy¬ 
domonas swimming in two and three dimensions [61]. 
Our case of colloids entrained at an oil interface of a sus¬ 
pension of swimming bacteria that are forming a FBI 
has several features that distinguish it from these stud¬ 
ies. First, the colloids, while in the incipient film, were 
in contact with the aqueous subphase and hence were 
coupled hydrodynamically to bacteria both in the film 
and in the bulk, a situation that is in some sense a hy¬ 
brid of the three-dimensional and quasi-two-dimensional 
systems considered previously. Second, as discussed be¬ 
low, FBI formation substantially impacts the rheology 
of the interface even at the earliest ages, with possi- 
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FIG. 1. (a) Ensemble-average mean-squared displacement of 
colloidal tracers during the active stage of FBI formation at 
an oil-water interface. For reference, the dashed line displays 
the mean-squared displacement of colloids at the oil-water 
interface in the absence of bacteria. The solid line displays the 
result of a fit to the data with Eq. 1. (b) Time autocorrelation 
function of the direction of instantaneous tracer displacements 
during the active stage of film formation. The line shows 
the result of an exponential fit to the data at large times 
{t > 0.03 s). 


ble qualitative consequences for the coupling between 
the swimming bacteria and the colloidal tracers in the 
film. Finally, Chlamydomonas are “pullers” while Pseu¬ 
domonas are “pushers”, and this distinction has conse¬ 
quences for the collective behavior and resulting hydro¬ 
dynamics of suspensions. Given these distinctions, the 
strong correspondence between the PDFs in Fig. 2(b) 
and those from bulk Chlamydomonas suspensions sug¬ 
gests that such self-similar distributions with Gaussian 
and Laplacian components might emerge more generi- 
cally in non-equilibrium active systems than previously 
thought. 

Another feature of our study of the active stage was the 
density of colloidal probes at the interface, which was 
large enough that correlated motion among the colloids 
could be observed, thereby providing information about 
the spatial correlations of the non-Brownian “kicks” the 
swimming bacteria impose. As an illustration, the inset 
of Fig. 3 depicts n for the colloids in the microscope field 
of view at one instant during the active stage. Align- 
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FIG. 2. (a) Normalized probability distribution functions for 
colloidal displacements at three lag times during the active 
stage of FBI formation. The dotted lines display the results 
of fitting Gaussian distributions in the regions of the peaks 
of the distributions, highlighting the enhanced, non-Gaussian 
probability for large displacements, (b) Normalized probabil¬ 
ity distribution functions from (a) plotted against normalized 
displacement (scaled with the root mean-squared displace¬ 
ment at each lag time), illustrating the collapse of the distri¬ 
butions onto a universal curve. The dashed line displays the 
result of fitting a Gaussian distribution in the regions of the 
peak. The solid line displays the result of a £t using the form 
given by Eq. 4. 

merit between the direction of motion of nearby colloids 
is clearly apparent. This coordinated motion is quanti¬ 
fied in Fig. 3, which shows the normalized pair direction- 
direction correlation function, 

Cnir) = {n,{T) (5) 

where the brackets represent an average over all pairs of 
particles i, j separated by distance r. The spatial cor¬ 
relations in tracer motion decay exponentially with sep¬ 
aration, as depicted by the solid red line in Fig. 3, which 
is the result of an exponential fit to the data with corre¬ 
lation length of 4.7 ± 1.7fim. Together, G„(t) and G„(r) 
give a quantitative picture of the short-time, spatiotem- 
poral correlations in tracer motion that characterize the 
dynamical behaviour of the interface during the initial, 
active stage of FBI development. 

Viscoelastic Transition 

Typically, the initial active stage of the film persisted for 
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FIG. 3. Normalized pair-direction correlation function as a 
function of the distance between colloids during the active 
stage of FBI development. Inset: Map of a section of the 
interface during the active stage of film formation showing 
the direction of motion of the colloidal tracers at an instant in 
time. Each colloid is represented by an arrowhead indicating 
its instantaneous direction of motion. 


less than 5 minutes, after which no motile bacteria were 
observed either at the interface or in the near-interface 
bulk. We ascribe the limited duration of the bacteria 
motility to the lack of nutrient in the suspension. The end 
to the active stage was reflected in a qualitative change to 
the probe dynamics in which the probe mean-squared dis¬ 
placement changed from superdiffusive at short lag times 
to subdiffusive. Once the bacteria ceased to move visibly, 
we treated the interface as a passive system close to ther¬ 
modynamic equilibrium and considered the probes to be 
undergoing thermally-driven Brownian trajectories from 
which the film rheology could be inferred. 

As mentioned above, the active stage was not always ob¬ 
served. Significantly, the ensuing mechanical changes of 
the interface, as inferred from probe mobility, appeared 
qualitatively independent of whether it was preceded by 
an active stage. Furthermore, as discussed below, the 
evolution in film rheology persisted for many minutes af¬ 
ter the end of the active stage, suggesting that the pres¬ 
ence of the active stage had limited impact on subsequent 
FBI evolution. From these observations we conclude 
that the film formation was primarily the consequence 
of polysaccharides and surface-active moieties produced 
by the non-motile (resting) bacteria. However, we can¬ 
not discount the possibility of some subtle effects to film 
formation due to biomixing by the swimming bacteria in 
those instances with an active stage. 

Here, we present results for the viscoelastic evolution of 
the FBI from an experiment in which no protracted ac¬ 
tive stage was visible at the start of film formation. We 
focus on this data set because (i) the population of (non- 
motile) bacteria at the interface was limited, leaving a 
relatively unobstructed interface for the colloids and fa- 
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cilitating the analysis of their Brownian motion, and (ii) 
the colloidal motion was less likely to be affected by any 
residual activity than in a trial with a fully realized ac¬ 
tive stage. We emphasize again that the trends observed 
were qualitatively the same as those in trials in which the 
viscoelastic film development was preceded by an active 
stage. 

Fig. 4(a) shows the ensemble-averaged mean-squared dis¬ 
placement (Ar^(t)) of the colloidal probes at several ages 
ta of FBI formation following creation of the interface. 
Because each mean-squared displacement is determined 
from 1 minute of video, the maximum lag time t is re¬ 
stricted to t < 2s to assure adequate statistics. Again 
for reference, the mean-squared displacement of colloids 
diffusing at the oil interface of pure Instant Ocean con¬ 
taining no bacteria is also shown. At our earliest mea¬ 
surement of film formation, the colloidal mean-squared 
displacement varies sublinearly with lag time, indicating 
a drag due to the development of a film with viscoelastic 
character. 

Within this limited dynamic range of accessible lag times, 
(Ar^(t)) is approximated as a power law, (Ar^(t)) ~ t” 
, with n < 1. As shown in Fig. 4(b), the power-law 
exponent n decreases steadily with increasing age, signi¬ 
fying increasingly subdiffusive motion. While in princi¬ 
ple the probe mobility is affected by both the interfacial 
film and drag from the surrounding bulk oil and water, 
given the large difference between (Ar^(t)) in the pres¬ 
ence or absence of the forming FBI, we can safely infer 
that the bulk contributions to the drag are insignificant. 
In this case, under appropriate conditions one can ob¬ 
tain the frequency-dependent interfacial shear modulus, 
G*{lS) = G'{u}) + iG"{oj), from the Brownian motion of 
the probes through a two-dimensional version of a gen¬ 
eralized Stokes-Einstein relation [62, 63], 

G*{uj) = - ,,, , (6) 

ttiujTu [(Ar2(t))] 

where Tu [(Ar^(<))] is the unilateral Fourier transform 
of the mean-squared displacement. Following Eq. (6), 
power-law behaviour in the mean-squared displacements 
of the colloidal probes, ^Ar^(t)) ^ t”, with n < 1, im¬ 
plies the film’s shear modulus has power-law frequency 
dependence G'(uj) ~ G"{uj) ^ w" [63]. Such weak 
power-law frequency dependence of G* is a character¬ 
istic of the rheology of a broad range of disordered 
complex fluids including concentrated microgel solutions 
[64], foams [65], paint [66], intracellular matrix [67], com¬ 
pressed emulsions [68], clay suspensions [69], and liquid- 
crystal nanocomposites [70], and is indicative of a broad 
spectrum of relaxation times. In most cases, the power- 
law exponent typically lies in the range n ^ 0.1 to 0.3. 
The soft glassy rheology model explains this response as 
a general consequence of structural disorder and metasta¬ 
bility, and provides a unifying theoretical framework for 
this behaviour. In this model, n serves as an effec¬ 
tive noise temperature, with systems approaching a glass 
transition as n —>■ 0. Thus, the steady decrease in n 



FIG. 4. (a) Ensemble-average mean-squared displacement of 
colloids at several ages during FBI formation at an oil-water 
interface for a trial with no active stage. For reference, the 
dashed line depicts the mean-squared displacement of colloids 
at the bare oil-water interface in the absence of bacteria. The 
solid lines display the result of power-law fits, (b) Power-law 
exponent characterizing the ensemble average mean-squared 
displacements, (Ar^(t)) ~ t" , of colloids at the oil-bacteria 
solution interface as a function of the age since formation of 
the interface. Inset: The interfacial elastic shear modulus Go 
at late ages, where the interface behaves like an elastic film. 
The elastic modulus grows logarithmically with age. 


with layer age reported in Fig. 4(b) points to increasingly 
glassy dynamics characterizing the structural response of 
the film. 

An important property of soft glassy systems is their 
non-equilibrium behaviour and spatial heterogeneity. As 
a measure of these features. Fig. 5(a)-(c) show the PDFs 
of colloidal displacements at three ages during the vis¬ 
coelastic transition. In each case, the PDF is shown 
at three lag times normalized by the mean-squared dis¬ 
placement at that lag time. At the earlier two ages, 
{ta = 57 and 2160 s), during which the viscoelastic char¬ 
acter of the film is evolving rapidly, the PDFs show a 
pronounced non-Gaussian component corresponding to 
enhanced probability of large displacements. These non- 
Gaussian contributions resemble those characterizing the 
colloidal dynamics in the active stage (Fig. 2); however, 
their origin in this case is different. Unlike at the active 
interface, the perturbations here are thermal, and each 











individual particle’s displacements are Gaussian. The 
non-Gaussian distributions result from variation in the 
mobility of different particles, evincing a spatially het¬ 
erogeneous interface of rheological microenvironments. 
Surprisingly, this heterogeneity diminishes at late ages, 
as illustrated by the closer-to-Gaussian distributions in 
Fig. 5(c), when the interface’s evolution has slowed and 
the film is nearly elastic. An interesting future study 
would be to compare this spatial heterogeneity with that 
of films formed in the presence of an extended stage of 
activity by swimming bacteria to investigate the role of 
biomixing in suppressing such heterogeneity. 

At late ages, when the layer behaves like an elastic film, 
and (Ar^) asymptotes to a constant value, we can extract 
an interfacial elastic shear modulus Gg from 


G' = 


keT 


{Ar^{t — )• oo)) 


(7) 


As shown in Fig. 4(b) inset, Gg grows logarithmically 
with layer age. Such logarithmic growth of the elastic 
modulus is a characteristic of the aging behavior of disor¬ 
dered soft solids such as colloidal gels and pastes [71, 72], 
indicating again the similarity between the rheology of 
the biofilm and that of soft glassy materials. 



B. Pendant Drop Elastometry 

A pendant drop of hexadecane is placed in contact with 
a bacterial suspension. As bacteria adhere to the drop 
interface, the FBI forms and the interfacial tension de¬ 
creases (Fig. 6). In early stages, the colloidal probe data 
above indicate that the interface is fluid, and the decrease 
in interfacial tension is attributable to surface active 
polysaccharides secreted by adsorbed bacteria. Eventu¬ 
ally, a film forms at the interface with significant elastic¬ 
ity. To probe this elasticity, the biofilm is deformed from 
its unstrained state by injecting or withdrawing liquid 
from the droplet, thereby changing the drop volume and 
area. Here we analyse a data set which captures key fea¬ 
tures of the FBI in late stages of development. A drop of 
hexadecane aged in a bacteria suspension for 22 hours has 
an apparently smooth surface with a shape corresponding 
to that for a pendant drop subject to an isotropic tension 
of 7 o = 21.0 ^^jm- When the drop volume is decreased, 
the drop shape changes dramatically, as shown in Fig. 7. 
For small enough compressions < 0.3), the drop 

shapes are used to infer the two-dimensional elastic mod¬ 
ulus K 2 D and the two-dimensional Poisson ratio V 2 d of 
the FBI by assuming a 2D Hookean constitutive model. 
For larger compressions, the drops wrinkle at their necks; 
this wrinkling phenomenon is discussed below. For still 
larger compressions, the drop resembles an empty, wet 
plastic sack with pronounced wrinkles. To analyse the 
small deformation regime, we integrate the anisotropic 
Young-Laplace equation for 2D Hookean films with un¬ 
strained tension corresponding to the apparent tension 


FIG. 5. Normalized probability distribution functions for col¬ 
loidal displacements at three lag times and three film ages 
(a) 57 s, (b) 2160 s, and (c) 88920 s during the viscoelastic 
stage of FBI formation. The solid lines display the results 
of fitting Gaussian distributions in the regions of the peaks 
of the distributions, highlighting the enhanced, non-Gaussian 
probability for large displacements. 



FIG. 6. Dynamic surface tension of a hexadecane drop aged 
in a bacteria suspension. Inset image shows the drop shape 
at the end of the aging cycle. The needle is 0.9081 mm wide 
and serves as a scale bar. 
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of the drop prior to compression, 70 = 20.44 jm- Con¬ 
sider the tension, strain and shape profiles, shown in 
Fig. 8 ; in this figure, lengths are scaled with the radius 
of curvature of the unstrained drop i?o = 1-30 mm and 
tensions and are scaled with 70 . The shapes of the 
unstrained (dashed curve) and compressed drops (solid 
curve) are shown in Fig. 8 (a). The apex of the com¬ 
pressed drop is located at the origin; this contour ends 
where the drop would intersect the needle at the right. 
The meridional tension Ts and parallel tension profiles 
are shown in Fig. 8 (b), for ~^^Iao = 0.13 for material 
parameters that capture the experimental drop profiles 
well, i.e. K 2 D = 21.3 V 2 D= 0.78. For much of the 

contour, the tensions are equal but non-uniform; they 
are highest at the apex where strains are weakest, and 
decrease steeply where the drop neck meets the needle. 
The high strain in this region is consistent with the area 
dilatation profile (Fig. 8 (c)). This anisotropic strain envi¬ 
ronment is key to the formation of wrinkles in the layer at 
large compressions [47]. In order for the film elasticity to 
give pronounced changes in drop shape like those evident 
in our experiments, K 2 D must be similar in magnitude to 
7 o. Values for this elastic modulus are obtained by mini¬ 
mizing an objection function, which is the cumulative dis¬ 
tance between edge points on the experimental drop and 
their respective nearest neighbor on the numerical drop, 
with respect to Bond number. Bo, location of the axis 
of symmetry, radius of curvature at the apex, Rq, elas¬ 
tic modulus K 2 d, and Poisson ratio, V 2 d- Examples of 
an intermediate fit and a fully optimized numerical drop 
contour are shown in Fig. 9. The resulting data, (Fig. 10), 
show K 2 D of a characteristic magnitude, in the range of 
10—27 with a downward trend suggesting that this 

is not a Hookean film in this strain regime. This modulus 
is far larger than that found for hexadecane drops emul¬ 
sified and aged over minutes in the presence of bacteria 
characterized by conventional pendant tensiometry and 
micropipette aspiration [33, 34], for which a modulus of 
5 '^^/m is reported. Pendant drop elastometry assuming 
isotropic films was performed for bacteria grown on min¬ 
eral and medium chain triglyceride oils, with elasticities 
ranging from 10 — 30 and 2—10 respectively 
[31]. Interestingly, the modulus we find is comparable to 
the spring constants for individual planktonic gram nega¬ 
tive bacteria as characterized by AFM, which range from 
10 — 20 ™^/m; bacteria may form biofilms with elastic 
moduli comparable to their own elasticities [73]. 

Note that values for the 2D Poisson ratio extracted by the 
scheme are highly variable (inset in Fig. 10). This vari¬ 
ation is not physically significant; we find that pendant 
drop shapes are very weakly dependent on this quantity. 
To illustrate this concept, a family of simulated profiles 
of Hookean pendant drops is shown (Fig. 11), generated 
from a single unstrained state and elastic modulus but 
differing 2D Poisson ratios. The profiles superpose for 
most of the drop contour, differing only in the high strain 
region of the drop neck. This is also the site of highest 
strain anisotropy where the role of the 2D Poisson ratio 


would be most significant. Well-resolved edge detection 
in the region of the neck would be extremely important 
to distinguish between these values; in many cases, in ex¬ 
periment, a droplet might attach to a needle before these 
curves peel away from each other, limiting the degree of 
precision with which one can evaluate the Poisson ratio. 
Furthermore, in our study, bacteria aggregates adhered 
to the FBI obscure the drop shape near the neck, pre¬ 
cluding the relevant analysis, which rely on having a well 
resolved drop profile. 

For large compressions, the neck region of the drop even¬ 
tually roughens. Upon continued compression, the rough 
regions deepen into wrinkles which lengthen vertically 
near the neck. Such wrinkled elastic films formed by FBI 
have been reported previously [27, 34]. 

The appearance of these wrinkles can be used to show 
that the 2D Poisson ratio of the film is positive. Fur¬ 
thermore, we can use these wrinkles to show the bending 
energy Eb of the film. Under strong compression, ten¬ 
sion reduces more in the parallel direction than in the 
meridional direction, which must support gravity’s pull 
in the vertical direction. These anisotropic tensions and 
associated strains cause wrinkles to form aligned with 
the meridional direction. The onset of this deformation 
occurs when the energy for compression becomes com¬ 
parable to the energy required to wrinkle, characterized 
by the bending energy Eb- At onset, wrinkles of a sin¬ 
gle wavelength Ac and associated amplitude should ap¬ 
pear for Ac ^ i?o which is related to the bending en- 

A4 

ergy according to Eb ~ ie^ 2£2 where is the length 
of wrinkled region and Tg is the meridional tension at 
onset [47, 74]. In our experiments, when wrinkles be¬ 
come apparent, however, deep wrinkles with wide spacing 
coexist with wrinkles of smaller amplitude and shorter 
wavelengths. The deep, widely spaced modes appear 
more prominent as compression is increased (Fig. 7(b)). 
Strictly speaking, the emergence of this band of wave¬ 
lengths indicates compression beyond the onset condi¬ 
tion. However, to infer an order of magnitude for the 
bending energy, we associate the longest wavelength with 
Ac (given as a range because of the decrease in drop ra¬ 
dius along the arc length), and estimate Tg from the last 
drop silhouette before wrinkling, to find an estimated 
range for As ~ 2 — 50 x 10“^® N ■ m [47]. This range 
of values is consistent with those reported previously for 
cross linked films of polymers (OTS) at fluid interfaces 
[47, 75, 76]. 

Finally, we investigate the hysteresis of these films under 
compression and re-expansion. The FBI stiffen signifi¬ 
cantly even after a single compression/re-expansion cy¬ 
cle. Fig. 12 shows this phenomenon. Panel A corresponds 
to the film-covered drop prior in its unstrained state. In 
this experiment, a pendant drop was aged in the bac¬ 
terial suspension for 67 hours, compressed to wrinkling, 
and subsequently re-expanded. Panel B is the drop after 
compression. Panel C corresponds to the drop after re¬ 
expansion to its initial volume and area, and D to areas 
well in excess of the initial drop (details are given in the 
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FIG. 7. Compression of FBI aged for 22 hours on a drop of hexadecane by withdrawal of oil. The compression increases from 
(a) to (c). The needle is 0.9081 mm wide and serves as a scale bar. As the drop is compressed, the film wrinkles. Approximate 
wrinkle analysis of (b), an intermediate stage of compression: Ac ~ 5 — 12 x 10“® m, L™ ~ 1.1 X 10“® m, r, ~ 5 


figure caption). Had the film not been hysteretic, the 
drop would have attained the same shape at a given vol¬ 
ume or area. In both panels C and the drop is appar¬ 
ently more spherical, illustrating its higher film tension 
in each, clear evidence of hysteresis. The microstructural 
origins of this apparent increase in tension have yet to 
be determined; possible mechanisms include irreversible 
folding of the film where wrinkles formed, effectively re¬ 
ducing the film area and creating stiffer regions, or irre¬ 
versible collapse of the porous microstructure of the film. 
These encapsulating films provide a strong barrier to co¬ 
alescence [33], in which a drop of oil aged in a bacte¬ 
rial suspension for 22 hours is compressed against a FBI 
formed at a planar oil-water interface. Significant com¬ 
pression is supported by the films when pressed against 
each other, with no evidence of rupture. 


IV. CONCLUSIONS 

We study the formation of films at interfaces of bac¬ 
terial suspensions (FBI), capturing associated changes in 
(micro)mechanics. Films may exhibit an early, active 
stage in which hydrodynamic interactions with bacteria 
near and at the interface create superdiffusive motion of 
colloidal probes in which the probe velocity has an ex¬ 
ponentially decaying memory. The probability distribu¬ 
tions of probe displacements are non-Gaussian but are 
self-similar at different lag times, even at times far ex¬ 
ceeding the characteristic time of the velocity memory. 
The active stage is followed by a viscoelastic stage dur¬ 
ing which the bacterial suspension forms a heterogeneous 
film with micromechanics that evolve over time. The film 
rheological properties behave similarly to a broad range 
of disordered complex materials in that the storage and 
loss moduli scale with w" where n < 1. Notably, the free 
interface displays remarkable changes in viscoelasticity 
after only a few minutes. This observation suggests that 
any study of the dynamics of bacteria suspensions with 



0.0 0.5 1.0 1.5 2.0 

z 


FIG. 8. Non-dimensional (a) theoretical shape, where (0,0) is 
the apex of the drop, and the needle would be at the ter¬ 
mination on the right, (b) principal tensions, parallel, r^, 
and meridional, Ts, and (c) local area dilatational (5A(s)= 
— 1) profiles of a drop under compression: Bo = 0.18, 
70 = 20.44 K2d = 17.0 = 0.78, Ra = 

1.24 mm, V = 9.8 fj.1. 


free interfaces must be treated with care, as collective 
dynamics depend sensitively on the mobility of bounding 
surfaces. Over time, the FBI become elastic via the for¬ 
mation of heterogeneous patches which eventually cover 
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FIG. 9. Fitting of experimental drop shape to numerical 
anisotropic Young-Laplace equation solution, (a) Experimen¬ 
tal contour and the numerically computed drop shape for an 
initial guess for the elastic modulus, K 2 D, and Poisson ratio, 
1^20 ■ (b) Best fit numerical solution to experimental data. 
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FIG. 10. Anisotropic analysis of a compression cycle: elastic 
modulus over 60 bins, with an average K 2 D = 20.81 
Inset: 2D Poisson ratio varies across the whole range of pos¬ 
sible values with no discernible trend. 


the interface. In this state, the film exhibits signature 
behaviour of thin elastic shells with associated bending 
energies and anisotropic tensions. We use pendant drop 
elastometry to obtain the properties of these films af¬ 
ter more than 20 hours of interface aging. By weakly 
compressing and expanding a FBI-covered oil drop, we 
determine the elastic moduli, which are comparable to 
’’spring constants” determined for Gram negative bacte¬ 
ria by AFM. By strongly compressing the FBI-covered 
oil drop, the films wrinkle, implying a positive Poisson 
ratio associated with the interfacial film. From the wave 
lengths of the wrinkles, we also extracted a bending en¬ 
ergy consistent with literature values for thin polymer 
films. Additionally, these FBI exhibit significant hys- 


.... i/= —0.9 12 = 0.9 12 =—0.36 



r [mm] 

FIG. 11. Illustration of weak dependence of pendant drop 
shapes on the Poisson ratio in the region of the drop neck on 
a drop with Bo — 0.18, 70 = 20.44 K 2 D ~ 17.19 

Ra = 1.22 mm. The needle is 0.9081 mm wide and serves as 
a scale bar. Inset plot shows an expanded plot of the family 
of V curves at the neck. Inset image shows aggregation by 
bacteria at the neck of the experimental drop. 

teresis upon expansion and compression. 

Finally, we have studied the behavior of a particular 
bacterium. Pseudomonas sp. ATCC 27259, strain P62. 
Different species or strains have different propensities to 
form FBI; these propensities depend on the presence or 
absence of nutrient sources. Furthermore, for a given 
strain, film (micro)mechanics likely depend on the sur¬ 
face structures and motility of the bacteria. We will ex¬ 
pand this work to address such issues, and to investigate 
community dynamics of competing and synergistic bac¬ 
teria. 
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FIG. 12. Compression and expansion of FBI on a drop of hexadecane, analysis performed with isotropic Young-Laplace 
equation; (a) initial drop, 70 = 12.7 V = 9.3 fil, (b) drop after full compression, (c) inflation to initial drop size, more 

spherical, which is indicative of higher tension, 70 = 21.5 ^^/m, V — 9.3 i-il, (d) increase in drop size, increased tension, 
70 = 29 V = 12.8 fil. The needle is 0.9081 mm wide and serves as a scale bar. 
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